Temperature evolution of physical parameters in the Inert 

Doublet Model 



Dorota Sokolowska 
University of Warsaw, Faculty of Physics, Warsaw, Poland 



January 19, 2013 



Abstract 



Inert Doublet Model is a minimal extension of the Standard Model with the second scalar 
doublet that may provide a Dark Matter candidate. In this paper we consider different variants 
of the evolution of the Universe after inflation, that lead towards the Inert phase today. We 
extend our previous analysis, in particular by discussing the co-existence of minima and pro- 
viding numerical examples of the evolutions of the mass parameters for the different types of 
the history of the Universe. We take into account the existing constraints, including the relict 
density data. 

1 Introduction 

The Inert Doublet Model (IDM) [Ij is a Z2 symmetric 2HDM, which for a special set of parameters 
may provide the Dark Matter (DM) candidate. The model contains two scalar SU{2) doublets. 
One is a "standard" scalar (Higgs) doublet $5, responsible for electroweak symmetry breaking and 
masses of fermions and gauge bosons as in the Standard Model (SM), with a Higgs particle hs- The 
second one is a "dark" scalar doublet which does not receive vacuum expectation value (v.e.v.) 
and does not couple to fermion^ In the model the following discrete symmetry of the Z2 type, 
which we call the D -symmetry, is presentFj 



All the components of the dark doublet are realized as the massive D-scalars: two charged 
and two neutral ones Dh and Da- By construction, they possess a conserved multiplicative 
quantum number, the D-parity, and therefore the lightest particle among them can be considered 
as a candidate for the DM particle. 

In this paper we discuss the evolution of the Universe during its cooling down after inflation, 
following the approach presented in [21 El HI E]. As before, we assume that the current state of 
the Universe is described by IDM. In this analysis we include all existing constraints on the model, 
together with the corresponding energy relict density, calculated by us [6J using micrOMEGAs j7]. 

In sec. 2 we list the basic properties of IDM and discuss the relevant astrophysical and collider 
constraints. Sec. 3 contains the summary of the basic assumptions of our approach and extended 
analysis of the possible today's states and types of evolution, as compared to [S]. In particular we 
trace the co-existence of minima. In sec. 4 we provide numerical examples of temperature evolution 
of physical parameters, like masses. 

^Our notations are similar to those in the general 2HDM with the change $i — > $5, $2 — > 
"^SM denotes SM fermions tpf and gauge bosons. 
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2 IDM 



2.1 Lagrangian 

We consider an electroweak symmetry breaking (EWSB) via the Brout-Englert-Higgs-Kibble (BEHK) 
mechanism described by the Lagrangian 



H 



T-V. 



(2) 



Here, describes the SU{2) x f/(l) Standard Model interaction of gauge bosons and fermions, 
which is independent on the reahzation of the BEHK mechanism. In the considered case the scalar 
Lagrangian Ch contains the standard kinetic term T and the D-symmetric potential V, which can 
describe IDM, with two scalar doublets $5 and 



V = -- 



+ A3 (*D*d) + A4 (^d'^s) + ^Ag V (*D$s)" 



(3) 



All parameters in V are taken to be real, with A5 < [5]. For the further discussion, it is useful to 
introduce the {fii, (12) phase space with: 



(4) 



Positivity conditions imposed on the potential guarantee the existence of the stable vacuum. 
They assure that the potential is bounded from below, meaning that the extremum with the lowest 
energy will be the global minimum of the potential (vacuum). The positivity constrains relevant for 
this analysis are: 



R 



Ai>0, A2 > 0, i?+l>0, 
A345/ A/A1A2, A345 = A3 + A4 + A5. 



(5) 
(6) 



Cy describes the Yukawa interaction of SM fermions ipj with only one scalar doublet $s, having 
the same form as in the SM with the change $ — > $5 (Model I for a general 2HDM). For quarks it 
reads: ^ 

Cy = QiT^sdR + Ql^^sUr + {h.c), (7) 

where is doublet of left-handed quarks, d^^UR are the right-handed quarks, F and A are 3x3 
matrices of Yukawa couplings in generations space, $ = ia2^* ■ Similar Yukawa term is introduced 
for leptons, namely L^T^s^R + {h.c). jCy respects D-symmetry in any order of the perturbation 
theory. 



2.2 Inert vacuum state 

Inert extremum, denoted by Ji, is realized if the extremum conditions for the potential V respect 

v"" = fii/^/X~i. (8) 
Ii realizes vacuum if following conditions are satisfied |5]: 

yUi > for any R, Hi > ^2 for R> 1, R^i > ^2 for < 1- (9) 

There are four dark scalar particles Dh, Da, and the Higgs particle hs which interacts with 
the fermions and gauge bosons just as the Higgs boson in the SM. 
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Inert state is invariant under the ^-transformation just as the whole basic Lagrangian g does. 
Therefore, the D-parity is conserved and due to this fact the hghtest D-odd particle is stable, being 
a good DM candidate. 

Masses of the scalar particles are: 

Ml = X,v^ = ml, , Ml, = h^^^l^ , 

A A A +A 

Mh^ = Ml, + , Ml^ = Ml, + . 

Assuming, as usual, that DM particles are neutral, we consider such variant of IDM, in which 

Mb±, Md^ > Md„. (11) 

After EWSB parameters of the potential V can be expressed by the four scalar masses and two 
self-couplings A345, A2 between the neutral scalar particles. Triple and quartic couplings between 
SM-like Higgs hs and DM candidate Dh-, i.e. DuDuhs and DnDnhshs, are proportional to A345. 
The second coupling, A2, is related only to a quartic self-coupling, DhDhDhDh. The remaining 
self-coupling, A3, governs the charged scalars' interactions: D^D^hs and D^D^hshs- 



2.3 Constraints 

Various theoretical and experimental constraints apply for the IDM (see e.g. |8]-|12j). 

The value of A345 strongly affects the DM interactions relevant for the DM energy relict density 
^Duh^- In general, for larger IA345I the relict density decreases due to the enhanced DhDh anni- 
hilation via s-channel Higgs exchange. This parameter also plays an important role in the indirect 
detection of DM [9l [10] and larger IA345I gives an enhanced flux of neutrinos and gamma rays. 

The value of the remaining coupling, A2, does not influence the DM relict density explicitly and 
so this parameter is usually fixed to arbitrary small value during the DM analysis of IDM [HI ITU] . 
However, value of A2 limits the value of A345 via the positivity constraints (A2 > IA345I/1/A1) and 
therefore it plays an important role in the analysis, as discussed in [6]. 



Collider constraints on scalars' masses Electroweak precision tests constrain strongly physics 
beyond SM. For IDM both light and heavy Higgs particle is allowed [Ij. Constraints on the mass 
splittings 6a = ~ ^Dhi ^± = ^d± ~ have been discussed in |9J. For a light Higgs boson, 
the allowed region corresponds to 5± ~ 5a with mass splittings that could be large. For heavy SM 
Higgs large 5± is needed, while 5a could be small. In this work we limit ourselves to the light SM-like 
Higgs boson. 

As D^,Da and Dh do not couple to fermions, the LEP limits based on Yukawa interaction for 
the standard 2HDM don't apply. However, the signatures are similar to neutralinos and charginos 
interactions in MSSM and the absence of a signal at LEP II was interpreted within the IDM in paper 
|llj . This analysis excludes the following region of masses: Mjj^ < 80 GeV, M^)^ < 100 GeV and 
5a > 8 GeV. For 5a < 8 GeV the LEP I limit Md^ + Md^ > Mz applies. 



Constraints on self-couplings The positivity constraints are imposed directly on quartic pa- 
rameters in the potential 
A cannot be large. The bound (called perturbativity constraint) is set typically to 

|A|<47r. (12) 

Astrophysical estimations of the energy relict density QoAih^ may be used to give the limitations 
for IA345I depending on the chosen value of masses of Dh and other scalars [9l [10]. However, they 
do not directly constrain the remaining quartic coupling A2. 



5||6). If we want to assure the perturbativity of the theory, self-couplings 
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DM relict density constraints The DM energy density in the Universe is estimated to [13j : 



(]^^^/i2 = 0.112 ±0.009. (13) 

In this analysis we assume that is a dominant component of the observed DM and density 
(13) is today's density of D^- 

Various studies [9^, ITI]] and our independent analysis show that for IDM there are three allowed 
regions of Md'- (i) light DM particles with mass close to and below 10 GeV, (ii) medium mass 
regime of 40 — 80 GeV and (iii) heavy DM of mass larger than 500 GeV. For purpose of this paper 
we concentrate on the medium mass region with the chosen set of masses (sec. 4). With this choice 
parameters A2 and A345 are free (up to the limitations discussed above). We consider different 
variants of their choice corresponding to different types of the evolution of the Universe. 

For our analysis use of standard available tools, ie. micrOMEGAs, is in fact limited, as this 
program neglects temperature dependence of physical parameters and a possibility of more than 
one phase transition. In paper |5] we concluded that if in the past there were sequences of phase 
transitions, then the Universe entered the inert phase with DM candidate at lower temperatures that 
in one-stage EWSB. This should be taken into account while solving the Boltzmann equations for 
DM relict density Ud- Furthermore, one should consider the latent heat of the first order transition 
or strong fluctuations in the sequences of the second order phase transitions |14| . For a single phase 
transition the use of present form of micrOMEGAs is more justified, however taking into account 
the evolution of masses and existence of different decay channels may provide corrections. In this 
sense, the relict density calculations in this paper should be considered as a preliminary estimate. 



3 Thermal evolution of the Universe 

We consider thermal evolution of the Lagrangian, following the approach presented in [U [31 E]. It 
allows to study the earlier history of the Universe after inflation. In the first nontrivial approximation 
the Yukawa couplings and the quartic coefficients A's remain unchanged, while the mass parameters 
m\ [i = 1,2) vary with temperature T as follow^ 

mUT)=ml-c,T', 

3A1 + 2A3 + A4 , 3g' + g" , g! + gl 3A2 + 2A3 + A4 , 3g' + g" ^^^^ 

Cl = \ \ z— ^, C2 = h 



6 8 2 ' ' 6 8 

Here g and g' are the EW gauge couplings, while gt and gt, are the SM Yukawa couplings for t 
and b quarks, respectiveljj^ 

In virtue of positivity conditions the sum of evolution coefficients is positive: C2 + Ci > 0. For 
R > both Ci > 0, while for R < arbitrary signs of Ci,2 are possible. In this work we limit ourselves 
to positive Ci,C2 as we consider only the restoration of EW symmetry for high T (this corresponds 
to the negative values of m^]^(T), m22{T) for high enough T) [15]. See |5] for more details. 

3.1 Possible minima during evolution 

As the Universe is cooling down the potential V (|3|), with temperature dependent quadratic coeffi- 



cients (14), may have different ground states. The general possible extrema are in form: 



i^s) =^^( a), i^D) =^=( ) , (15) 



with vs > and v'^ = vl + Iv^l + v?. 



s ^ I^d' ' 



•^Formulae for ci,C2 were obtained recently by G. Gil (Master Thesis, 2011); they correct the corresponding formulae 
given in |5]- 

"Normalization of couplings: 5, = V^TOj/w, {gt ~ 0.99,^6 « 0.02); g = 2Mw/v = 0.652, g' = 0.351. 
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Extrema and vacua 


name of extremum 


properties of vacuum 


vev's 


EW symmetric: EWs 


Massless fermions and bosons and massive 
scalar doublets. 


Vd = 0, Vs = 


inert: Ii 


Massive fermions and gauge bosons; scalar 
sector: SM-like Higgs hs and dark scalars 
Dh,Da,D^ with DM candidate Dh- 


Vd = 0, 
2 2 A'l 


inertlike: I2 


Massless fermions and massive gauge 
bosons; scalar sector: Higgs particle 
(no interaction with fermions), four 
scalars Sh, Sa, S"^, no DM candidate. 


Vs = 0, 

2 2 

Vh=v = n- 


mixed: M 


Massive fermions and bosons, 5 Higgs 
particles: CP-even h and H, CP-odd A 
and charged H^, no DM candidate. 


_ fil-Rfl2 ^ 

~ v^(l - m ^ 

_ fi2-Rfil ^ 



Table 1: General properties of the extrema and vacua, following [5j. 



Possible neutral solutions {u = 0) are: EW symmetric EWs, inert Ji, inertlike I2 and mixed M. 
Their general properties are summarized in table [T] [5j, see also the appendix A. 

There are three distinguish allowed regions of parameters X's, the best parametrized by the 
parameter R namely a) i? > 1, b) 1 > _R > 0, c) > i? > -1. 

The EW symmetric state [EWs) exists for every value of R if and only if both /ii < and 
IJ,2 < 0, being the only existing extremum (and thus the vacuum). For R> 1 (fig, la) the energy of 
the mixed extremum M (if it exists) is always higher than for the other extrema, so it cannot be 
the vacuum p]. Possible EWv vacua in this case are Ji and I2. Fig lb shows the allowed regions for 
< i? < 1. Again we have the regions of EWs, Ii and I2 vacua, but now also mixed extremum M 
can be realized as a vacuum. Case of i? < is similar to the previous one but with the wider region 
of M vacuum (figjlc|. Summary of existing extrema, local minima and vacua for various R regions 
can be found in table 2. 



3.2 Possible sequences of phase transitions 

We use {ni{T) , fi2{T)) phase diagrams and the redefined evolution coefficients Ci, C2 to determine 
the sequences of transitions between different vacua as T decreases: 

/ii(T) = mj,iT)/VX~,, /i2(T) = mUT)/VX~2, (16) 

Cl = Ci/a/A7, C2 = C2/a/A^, C = C2/C1. (17) 

Figs.([T^,b,c) show all possible types of evolution from the EW symmetric phase towards the inert 
phase today. The possible evolutions are represented by the rays directed towards the growth of 
time, ie. the decrease of the temperature. Today's values are defined by = fii{0),^2 = /U2(0) and 
marked by dots in figures. 

The relevant temperatures with the corresponding conditions for the phase transitions are |5]: 
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a) Extrema and vacua for R>1 


/i2 > and /i2 > yUi-R 


/2 is the vacuum. For /ii > Ji is an extremum, but not a minimum. 


< < ^2 < /Ul-R 


I2 is the vacuum (i.e. global minimum) and Ji - a local minimum. 


< fiiR^^ < /i2 < /Ui 


Ji is the vacuum (i.e. global minimum) and /2 - a local minimum. 


/^i > and fi2 < jiiR^^ 


Ji is the vacuum. For fi2 > I2 is an extremum, but not a minimum. 


b) Extrema and vacua for 1 > i? > 


/i2 > and /i2 > yUi-R~^ 


I2 is the vacuum. For fii > Ii is an extremum, but not a minimum. 


/ii > and /i2 < yUi-R 


Ji is the vacuum. For fi2 > I2 is an extremum, but not a minimum. 


< HiR < jj,2 < fiiR^^ 


M is the vacuum, Ji and /2 are the extrema. 


c) Extrema and vacua for > i? > — 1 


/i2 > and fi2 < l^iR'^ 


I2 is the vacuum, no other extrema. 


/ii > and /i2 < yUi-R 


Ji is the vacuum, no other extrema. 


yU2 > Max /^li?) 


M is the vacuum. Extrema: I2 for fi2 > 0, Ji for /xi > 0. 



Table 2: extrema and vacua for various R regions. 



Tews,i — \/ /^l/ Cl 
7£;ws,2 = -x/ /i2/ C2 



^2,1 



2,M 



A^l - /i2 
Cl - C2 

/il — R^2 
Cl - i?C2 ' 



{^x^TewsS) = 0), 

(/U2(7£;h/s,2) = 0), 



(/^l(^2,i; 

7m,i = 



{i^2{T2,m) = I^i{T2,m)/R, f^l{Tj 



/^2(T2,i)), 



-R/il — /i2 

Rci - C2 



(18) 
(19) 

(20) 
(21) 



M,l, 



fl2{TM,l)/R) 



3.2.1 R>1 

For i? > 1 the possible types of evolution that lead towards the inert phase today can be divided into 
two groups (fig, la, table [s]). The first one corresponds to a single phase transition of the 2nd-order 
at the temperature (18), shown by rays la, Ila and III. The second type of evolution is in form of a 



sequence of two phase transitions represented by the rays IV and V. For these rays, the first phase 



transition EWs — )■ I2 (EWSB) realized at temperature (19) is a 2nd-order phase transition, while 
the second one I2 — ?■ /i, realized at temperature (20), is a Ist-order transition. 



We notice that the considered R> 1 case contains an unique opportunities of co-existing minima 
(vacuum Ji and local metastable minimum I2) for rays III, IV and V. Also, only in this case there is 
a possibility of the Ist-order phase transition (rays IV and V). That cannot be realized for the other 
values of R. For ray IV the co-existence is temporary and the local minimum I2 disappear for low 
temperatures, while for rays III and V it still exists for T = 0. 

For R = 1 the phase space of co-existing minima Ji and I2 is reduced to line /ii(T) = fi2{T). 
Here rays III and V are not possible, unless additional condition of fii = ^2 is fulfiled. That however 
leads to the existence of two degenerate minima for T = 0. 
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(a) i? > 1 (b) 1 > i? > (c) > i? > -1 



Figure 1: Possible vacua and evolution to a current states (dots) represended by rays on (yUi,/i2) 
plane for a) i? > 1, b) 1 > i? > 0, c) > i? > —1. The boundary lines are A : fi2 = /ii-R, 
B : fi2 = fJ'i, C : fi2 = fii/R. Blue (dark shade) region represents I2 vacuum, yellow (light shade) 
region - Ji vacuum and green (medium shade) region - M vacuum. In the hatched regions between 
lines A, B and B, C Ii and I2 minima co-exist. 



If c 7^ 1 then after EWSB only one EWv minimum is created (either Ji or I2). However, there is 
a possibility that c = 1 and fii = fi2- In this case after EWSB two degenerate minima appear and 
they have the same energy for the whole time after symmetry breaking. Note, that because of the 
form of Yukawa interaction the fermionic contribution appears only in Ci and so c = 1 is possible 
only if Ai 7^ A2. 



Ray 


Sequence 


Conditions 


Co-existing min. and extr. for T = 


la 


EWs ^ h 


/i2 < 




Ila 


EWs -> h 


< /i2 < Min(/iic, /ii/2"^) 


I2 extremum 


III 


EWs h 


fiiR^^ < ji2 < Min (/iic, /ii) 


I2 local minimum 


IV 


EWs ^ h-^ h 


jl{c < jJ,2 < jliR^^ 


I2 extremum 


V 


EWs ^ h 


Max (/iic, /iii?^"*^) < /i2 < /ii 


I2 local minimum 



Table 3: Possible rays for i? > 1. 



3.2.2 0<i?<l 

In the < -R < 1 case for every point in phase diagram (/Xi(T), /i2(T)) there is only one existing 
minimum (and so it is a vacuum), as shown in fig, lb , and all transitions are of the 2nd-order. We can 
reach the today's inert phase by a single phase transition or through the sequence of three 2nd-order 
phase transitions (table |4j First type of sequence EWs — )■ Ji is realized by rays lb and lib, which are 
the analogs of rays la and lib. For ray VI, which corresponds to the sequence EWs — t- /2 — !■ M — )■ Ji, 
EWSB happens at temperature (19). Then there are two more transitions in this sequence: from I2 
into M and from M into the inert vacuum Ji, with the last transition at temperature (21). 
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Ray no. 


Sequence 


Conditions 


Co-existing min. and extr. for T = 


lb 


EWs h 


/i2 < 


- 


lib 


EWs -> h 


< /i2 < Min (/iic, fiiR) 


I2 extremum 


VI 


EWs h ^ M ^ h 


jjLiC < 1^2 < l^iR 


I2 extremum 



Table 4: Possible rays for 1 > i? > 0. 



3.2.3 



-1 < < 



We consider the restoration of EW symmetry in the past (both Ci, C2 positive) and there is only one 
possible ray Ic, which is similar to the rays la and lb (table Isl fig lc) [5J. 



Ray no. 


Sequence 


Conditions 


Co-existing min. and extr. for T = 


Ic 


EWs ^ h 


/i2 < /Uii? < 





Table 5: Possible rays for > i? > — 1. 



4 Temperature evolution of physical parameters 

In the previous sections to describe the history of the Universe we used the parameters of the 
Lagrangian, namely /ii,/i2 and i?, c. Relations between them gave us information about the possible 
sequences of the phase transitions. In this section we illustrate the underlying temperature evolutions 
of the physical parameters (i.e. masses of the particles) for various rays, each representing the 
different history of the Universe. Here, it is useful to fix six free parameters of the model in form of 
four physical masses 

Mhs, Mb^, Md„, Md± (22) 

and two self-couplings: A345 and A2. The considered values of masses and X's are chosen in agreement 
with existing constraints both from the colliders and the DM abundance measurements, as discussed 
in sec. 2.3, The values of Q^Mh"^ were calculated with the existing micrOMEGAs code. We expect 



the VtriMh to lie in the 3cr WMAP allowed range: 

0.085 < VtoKih'^ < 0.139. (23) 

As we treat those results as an estimate only the fact that for some rays the calculated value is 
slightly outside the this range does not exclude automatically those rays. 
The following scalar mass set was used for the today's inert phase: 

= 120 GeV, = 45 GeV, Md^ = 115 GeV, Md± = 125 GeV, (24) 

with self-couplings A345 and A2 different for each ray [18|^ 



Below we show mass evolutions arising from (14) represented by different rays for all three R 



regions: rays la. III, IV and V for i? > 1, ray VI for 1 > i? > and ray Ic for R < (figsj2H7j). We 
plot the temperature dependent masses of the scalars for every vacua that is realized for a chosen 
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The detailed discussion of the importance and constraints of those couplings is in progress [6J 
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ray. We also present the temperature dependence of the "mass parameters" for scalar states of local 
minima. In addition we present dependence on the temperature of two other physical parameters, 
namely the top quark mass mt{T) and v{T) (proportional to the W boson mass Mvi/(T)). 

The initial state of the Universe is the EW symmetric phase with two massive scalar doublets 
$5, and massless both fermions as well as gauge bosons. In this vacuum there are eight massive 
scalar states that come from the two scalar doublets, denoted on plots by ^s,^d, with masses equal 
to \mii(T)\/\/2, \m22iT)\/ \/2, respectively. After EWSB the Universe enters one of the symmetry 
violating vacuum with the proper particle content, with masses described by the formulas presented 
in the appendix |Xj 



4.1 The R> 1 case 
4.1.1 Ray la 



Figs. ([2p,,b) show mass evolution for the ray la, corresponding to the single phase transition EWs — )■ 
Ji, for the following set of parameters: 



Ao = 0.012, A 



345 



0.065 



ci = 0.913, C2 = 0.309. 



(25) 



The EWSB happens at T = 125.6 GeV and Universe enters the phase with massive fermions and 
gauge bosons (fig, 2b). The Universe stays in the Ii phase with the mass of DM candidate nearly 
constant Md„\t) ^ 45 GeV, fig. [2a} 




T[GeV] 




- T[GeV] 



(a) masses of scalar states in /i: M^^ (blue), M^)^ (red), (b) v{T) (solid) and mt{T) (dashed) for EWs (black) 

Mjj+ (purple), Mhs (green) and EWs: (black), <i>s and Ii (red) 

(grey) 

Figure 2: Evolution of masses (ray la): EWs — )■ Ji. 



4.1.2 Ray III 

Ray III (figs|3)i,b) can be realized for: 

As = 0.02, A345 = 0.117 ci = 0.93, C2 = 0.33. (26) 

Again here there is a single phase transition and after EWSB the Universe enters the /i phase 
at T = 124.5 GeV. This case is different from the previous one as at T = 57 GeV another minimum 



- local minimum I2 - appears (shaded region in fig, 3a). Dashed lines show the change of "mass 



parameters" for corresponding scalar states of this local minimum (see appendix [A| . 
4.1.3 Ray IV 

Figs.(|4^,b) show the mass evolution for the ray IV, for the following set of parameters: 

A2 = 0.068, A345 = 0.16 ci = 0.944, C2 = 0.368. (27) 
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The EWSB happens at T = 123.6 GeV when Universe enters the inerthke phase I2 with massless 



fermions and massive gauge bosons (fig, 4b). As the time grows another extremum appears, which 
later becomes a local minimum Ji. The first-order phase transition I2 — ?■ /i happens at T = 
123.1 GeV. Note, that two minima coexist during a period of time AT ^1.5 GeV (shown by the 



shaded region in fig, 4a). The "mass parameters" of the scalar states in the local minima /i and 
I2 are shown. The discontinuity in masses of physical particles: scalars, fermions and gauge boson 
(proportional to v) is visible. Universe enters the inert phase Ji with massive fermions, gauge bosons 
and scalars, among them with DM candidate Dh and their mass evolution continues up to the T = 
mass values. 

4.1.4 Ray V 

Figs.([5^,b) show ray V, which can be realized for the following parameters: 

A2 = 0.05, A345 = 0.17 ci = 0.948, C2 = 0.363. (28) 

First, there is EWSB into the I2 phase at T = 131 GeV. Then at T = 113.5 GeV the local 
minimum Ji appears. The first order I2 — h transition happens at T = 71 GeV, and I2 becomes a 
local minimum, which does not disappear up to T = 0. These two minima coexist during a period 



represented by the shaded region in figjSa} Dashed lines correspond to the "mass parameters" of the 
scalar states in the local minima Ji and 12- Note, that for this ray the final phase transition happens 
at the lower temperatures than in the other cases, where ratio Afo^^^ji/T is of the order 1. 



4.2 The 1> R>0 case 
4.2.1 Ray VI 



Figs. ([6p,,b) show the evolution of the Universe along the ray VI. At every temperature there is only 
one minimum, all three transitions are of the 2nd-order. We study mass evolutions for: 



A2 = 0.125, A, 



345 



0.17 



ci = 0.947, 02 = 0.40. 



(29) 



EWSB happens at T 



124.8 GeV. Universe enters the inertlike phase I2 with massless fermions 
121.1 GeV the mass of Sh particle goes to and the 



and massive gauge bosons (fig, 6b). At T 
2nd-order transition to the M phase takes place. This phase is very short lived, at the beginning and 
at the end the mass of h particle goes to zero, at T = 121.1 GeV and T = 120.9 GeV, respectively. 
At this last 2nd-order transition Universe enters the inert phase Ji with the DM candidate D^- The 
mass evolution continues to the today's values of masses. Note, that here rays V and VI have the 
same A345, while they differ by the value of A2. 













I 

- 






h 


EWs 












% 


Oh 


hs 





20 40 60 80 100 120 140 



T[GeV] 




(a) masses of scalars in /i : AIjj^ (blue) , Md^ (red) , Miy+ 
(purple), M/is (green) and EWs: (black), $5 (grey). 
Dashed lines - mass prameters of local 12- Shaded region - 
local minimum 12- 



T[GeV] 

20 40 60 80 100 120 140 

(b) v{T) (solid) and mt(r) (dashed) for EWs (black) 
and Ii (red) 



Figure 3: Evolution of masses (ray III): EWs — )■ /i. 
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4.3 The case > R > -1 



4.3.1 Ray Ic 

Ray Ic shown in figs.([7^,b) is realized for: 

A2 = 0.1, A345 = -0.115 ci = 0.852, ca = 0.293. (30) 

In this case there is a single phase transition (EWSB) at T = 130 GeV. Universe enters the inert 
phase Ji with massive fermions, gauge bosons and DM candidate Dh- 



4.4 Relict densities for considered rays 



We use the micrOMEGAs code [7J for a rough estimation of the relict density for the rays described 
in the previous sections 4.1-4.3. Results for VL^Mh^ are presented in figjs} where dots represent values 
obtained for the different rays. Ray la with VLoMh"^ = 0.31 is excluded by the WMAP results. For 
rays Ic and III values of Q£,Mh'^ are within the WMAP range. As IA345I grows the observed Q^Mh"^ 
decreases and for rays IV, V and VI it is below the lower WMAP limit. Q^Mh"^ for rays VI and V 
are equal, as those rays differ only by the value of A2, which does not enter explicitly the rates for 
processes relevant for the DM abundance. 

As discussed before, those values are an estimate of the real relict density. Sizable corrections 
are expected especially for ray V, because in this case the temperature of the final phase transition 
(Ist-order) is the lowest. 

Figjs] contains QoMh^ as a function of A345. Calculation was done for A2 = 0.15. Although 
floM does not depend explicitly of the value of A2, its value limits the range of A345 we can scan 
over because of the positivity constraints and necessary conditions for the inert minimum to be 
a global minimum [6]. In the considered case we can have the physical solutions only in range 
0.2 > A345 > -0.2. 

For the large mass splitting that we chose, the coannihilation is not important and the main 
process is DhDh bb. For low values of IA345I the relict density is high, as the annihilation via hs 
is low. It becomes more important and it lowers the resulting Q^mIi"^ as IA345I grows. We find that 
two A345 regions, A345 G (-0.105, -0.13), (0.105,0.13), are in agreement with the WMAP limit (pf. 
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T[GeV] 
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(a) masses of scalar states around phase transitions for Ii 
(thick Hnes) and I2 (thin Hnes). Green - {hs,h£)), purple 
- iD^,S^), red - Ida,Sa), blue - (Dh,Sh). Dashed and 
dotted lines - local minima parameters. Shaded region - 
co-existence of minima Ii and I2 ■ 
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(b) v{T) (soHd) and mt(r) (dashed) for EWs (black), 
I2 (blue) and Ii (red) 



Figure 4: Evolution of masses (ray IV): EWs I2 h- 
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(b) v{T) (soHd) and mt(r) (dashed) for EWs (black), 
I2 (blue) and Ii vacuum (red) 



Figure 5: Evolution of masses (ray V): EWs — )■ J2 — t- /i. 



5 Discussion and conclussions 



In this paper we studied further the evolution of the Universe after inflation towards the present 
inert phase. We considered different types of evolution, which may be parametrized by using the 
parameter R. 

We extend the studies with respect to the previous paper [5] by discussing the possibility of 
co-existence of the local minima. This oportunity may be realized for two types of rays - ray IV and 
V. In the first case the co-existence is only temporary, as the phase transitions happen in a short 
period of time and the local minimum I2 dissapears shortly after the Universe enters the inert phase. 

However, we find that also other situation is possible. For ray V the local minimum I2 does not 
disappear and it exists for T = 0. Also for ray III, which corresponds to a single phase transition, 
the local minimum I2 appears for the later stages of evolution and it exists for T = 0. Furthermore, 
this ray gives a good relict density. 

We also stress the fact that the intermediate phases I2, M for rays IV and VI are short-lived. In 
this sense those sequences may be considered as a similar to the one with the single phase transition, 
however the latent heat of the Ist-order phase transition may affect the evolution. 

One should keep in mind the limitiations of using micrOMEGAs in the context of our analysis. 
First, we consider only the two body final states and so, for our mass range, the actual Qdm may be 
lower [19]. Furthermore, the temperature dependence of masses and the latent heat from the first 
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(a) masses of scalar states around phase transitions for (b) w(T) (solid) and mt(T) (dashed) for £^Ws (black), 
Ii (thick lines), M (thin) and I2 (dashed). Green - I2 (blue), M (green) and Ii (red) 
(hs, H,hD), purple - (D±, iJ±, 5±), red - (Da, A,Sa), blue 
- {DH,h,SH) 

Figure 6: Evolution of masses (ray VI): EWs — t- /2 — ?■ M — t- Ji. 
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T[GeV] 



(a) masses of scalar states in Ii: M^,^ (blue), M^,^ (red), (b) v{T) (solid) and mt{T) (dashed) for EWs (black) 

Mu+ (purple), (green) and EWs: (black), $5 and Ii vacuum (red) 

(grey) 

Figure 7: Evolution of masses (ray Ic): EWs Ii. 




Figure 8: Relict density ^^dm^^('^345) (red curve) with WMAP limits (dashed lines). Dots correspond 
to particular rays. 

order phase transition will influence the value of VL£,Mh'^- We expect that the relict density for rays 
IV- VI will be higher, as the flnal phase transition happens at the lower temperature than the EWSB 
for the other rays. The corrections from the Ist-order phase transition will result in the splitting 
between rays V and VI, which for now have the same relict density. We should also take into account 
the fact, that during the evolution the relation between the mass of the DM particle M/j^ and masses 
of the possible decay states {M\y, Mt, M^) may vary. For certain rays decay channels other than hh 
may play important role and affect the energy relict density (for example, DhDh — > WW for 
Moh{T) > Mw{T)). The detailed analysis of those effects is in preparation |14] . 
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Ministry of Science and Higher Education Grant N N202 230337. 
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A Vacua properties 



The neutral solution of extremum conditions give the following values of VsjVo parameters: 



EWs 
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vd = 0, vs = 0, 
vd = 0, vl = v"^ = 
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11 . 



8Ai 



vs = 0, 



m: 



22 



S 



m. 



A2 ^ 8A2 

2 _ '^LAi ~ A345m^-^ 



22 . 



"^?iA2 - A345m22 



A1A2 — A^ 



345 



A1A2 



^345 



M 



m\^\2 - 2A345m^im22 + rn^^M 
8(AiA2 - AI45) 



(31) 
(32) 

(33) 
(34) 



Some of the equations (32)-(34) can give also negative values of f| or In such case the 
extremum, described by corresponding equations, is absent. 



A.l EW symmetric vacuum EWs 



The electroweak symmetric extremum EWs (31) is a minimum if mf^^ 22 < 0- Gauge bosons and 



fermions are massless, while the doublets have non-zero masses and 



I'"i2l 

V2 



, respectively. 



A. 2 Inert-like vacuum lo 



The inertlike vacuum I2 is "mirror-symmetric" to the inert vacuum /i, compare (32) and (33), with 
one Higgs particle h£, and four scalar particles: Sh, Sa, S^. The interaction among scalars and 
between scalars and gauge bosons are mirror-symmetric as well, so the only difference between I2 
and Ji arises from the Yukawa interaction. 

The inertlike vacuum I2 violates D-symmetry. The Higgs boson h couples to gauge bosons just 
as the Higgs boson of the SM, however it does not couple to fermions at the tree level. The scalars do 
interact with fermions. Therefore, here there are no candidates for dark matter particles. Note that 
all fermions, by definition interacting only with $5 with vanishing v.e.v. ($5) = 0, are massless. 
(Small mass can appear only as a loop effect.) 

The masses of the Higgs boson hn and S-scalars are given by (cf. ([lo))) with m^^^ ^ ^22- 



A. 3 Mixed vacuum M 

The mixed extremum]^ M violates the Z2 symmetry. In this vacuum we have massive fermions and 
no candidates for DM particle, like in the SM. There are five Higgs bosons - two charged and 
three neutral ones: the CP-even h and H and CP-odd A. Couplings of the physical Higgs bosons 
to fermions and gauge bosons have standard forms as for the 2HDM, with the Model I Yukawa 
interaction. 

Masses of scalars are as follows (see, e.g. [171 [2]) 

Mh. = -^^^v\ Ml = -v'K {v^ = vl + vl). (35) 
The neutral CP-even mass matrix is equal to 

>l=f/^^^ ^^ffA. (36) 

^Sometimes called a normal extremum N, see e.g. [16] 
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The mass matrix (36) gives masses of the neutral CP-even Higgs bosons: 



with sign + for the H and sign — for h. 
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